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ABSTRACT 


An analysis of the acoustical and flowfield environment for the scaled 1-pound-force (lbf) thrust 
tabletop motor was performed. This tabletop motor from NASA Stennis Space Center is com- 
posed of Plexiglas burning in gaseous oxygen with a graphite insert for the nozzle portion. The 
nozzle has a throat diameter of 0.2 inch and an exit diameter of 0.38 inch. With a chamber pres- 
sure at 55 pounds per square inch absolute (psia), a normal shock is formed immediately down- 
stream of the nozzle exit plane as the combustion products exhaust into the ambient at atmos- 
pheric pressure. The jet characterization is based on computational fluid dynamics (CFD) in con- 
junction with Kirchhoff surface integral formulation and compared with correlations developed 
for measured rocket noise and a pressure fluctuation scaling (PFS) method. Predictions and 
comparisons are made for the overall sound pressure levels (OASPL’s) and spectral dependence 
of sound pressure level (SPL). The overall objective of this effort is to develop methods for scal- 
ing the acoustic and flowfield environment of rockets with a wide range of thrust (1 lbf to 1 mil- 
lion lbf). 
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ACRONYMS, ABBREVIATIONS, AND SYMBOLS 


CFD 

computational fluid dynamics 

co 2 

carbon dioxide 

dB 

decibels 

FFT 

Fast Fourier Transform 

go 2 

gaseous oxygen 

Hz 

hertz 

kHz 

kilohertz 

KSC 

John F. Kennedy Space Center 

LaRC 

Langley Research Center 

Lbf 

pound force 

LNG 

liquified natural gas 

LST 

Launch Systems Testbed 

NASA 

National Aeronautics and Space Administration 

OASPL 

overall sound pressure level (ref. 20 p Pa) 

PFS 

pressure fluctuation scaling 

psia 

pounds per square inch absolute 

°R 

degree rankine 

RNP 

rocket noise program 

SPL 

sound pressure level (ref. 20 p Pa) 

Nomenclature 



c = sound speed, ^jyg c R u T IW , ft/s 
dj = nozzle exit diameter, ft 

F = total thrust, lbf 
/ = frequency, Hz 

g c = gravitational constant, (lbm-ft/lbf.s 2 ) 

/ = specific impulse (sec) 

M = Mach number, Vic 

p = static pressure, lbf/ft 2 

p' = acoustic pressure disturbance, p- p m , lbf/ft 2 

Re = Reynolds number, pu j d j / p 

R u = universal gas constant, lbf.ft/(lbm-mole.R) 

fj = nozzle exit radius, ft 

St = Strouhal number, fdj / u ] 

T = temperature, R 

u = axial velocity, ft/s 

u c = jet center-line velocity at exit, ft/s 
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ACRONYMS, ABBREVIATIONS, AND SYMBOLS (cont) 


V = velocity, m/s 

W = molecular weight, lbm/lbm-mole 
jc= axial distance from the nozzle exit, ft 
z = radial distance from jet axis, ft 


Greek Symbols 

ju = dynamic viscosity, lbm.ft/s 
p = density, lbm/ft 3 
Y = isentropic exponent 

O = fluctuating pressure spectral density, dimensionless 

O) = circular frequency, radians/s 

£2 = Strouhal number (pressure fluctuation method) 


Subscripts 

c = centerline 
e = nozzle exit 
j = jet 

t = stagnation or total, or nozzle throat 
CO = ambient 
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1. INTRODUCTION 


Acoustic loads in a launch vehicle environment represent a principal source of structural vibra- 
tion and may be critical to the proper functioning of vehicle components, ground support struc- 
tures, and equipment in the immediate vicinity of the launch pad. A knowledge of acoustic 
loads, including the overall sound pressure level (OASPL), sound pressure level (SPL) spectrum, 
and the distribution (or correlation) of surface acoustic loads is necessary to provide the input for 
vibroacoustic analysis and evaluation of structural integrity. In the design of launch vehicles, it is 
highly desirable that data on acoustic loads (near-field and far-field noise levels) be generated 
both analytically and from testing of small-scale and full-scale models. Since full-scale acoustic 
and vibration testing is often cost prohibitive, the option of small-scale testing combined with 
analysis methods remains as a practical alternative. Accurate characterization of acoustic loads 
on launch pad structures thus proves to be a formidable challenge. 

At NASA Kennedy Space Center, a Launch Systems Testbed (LST) is currently under develop- 
ment to establish a capability to simulate small-scale launch vehicle environment for use in test- 
ing and evaluation of launch pad designs for future space vehicles. The test program provides for 
the scale model testing of cold jets (nitrogen), hot inert jets (helium), and combusting jets 
(solid/liquid rocket fuels). Concurrently, prediction methods are also under development to com- 
plement the small-scale test program. These prediction methods include computational fluid 
dynamics (CFD), analytical correlations, and pressure fluctuation scaling (PFS) methods. 

Noise from subsonic jets is mainly due to turbulent mixing, comprising the contributions of 
large-scale and fine-scale structures (Lighthill, 1952, 1954). The turbulent mixing noise is 
mainly broadband. In perfectly expanded supersonic jets, the large-scale mixing noise manifests 
itself primarily as Mach wave radiation, caused by the supersonic convection of turbulent eddies 
with respect to the ambient fluid. In imperfectly expanded supersonic jets, additional noise is 
generated on account of broadband shock noise and screech tones. 

This report summarizes the analytical studies of a small-scale tabletop rocket that is being ac- 
quired from NASA Stennis Space Center for the purpose of establishing the initial test facility 
and instrumentation. This tabletop motor has Plexiglas fuel burning in gaseous oxygen (GO 2 ) 
with an off-design nozzle (overexpanded exit Mach number of 2.6 and 1 Ibf thrust). An impor- 
tant characteristic of the tabletop rocket is the existence of a normal shock in the immediate 
downstream of the nozzle exit, resulting in a subsonic jet structure. Comparisons of OASPL and 
SPL spectrum are presented for CFD simulations, rocket noise correlations, and the PFS method. 

2. DESCRIPTION OF THE TABLETOP ROCKET 

2.1 NOZZLE GEOMETRY 

Figure 1 shows a schematic of the nozzle geometry for the tabletop motor, as obtained from 
NASA Stennis Space Center (Field 2001). Both the convergent and the divergent portions of the 
nozzle are conical. The nozzle has a throat diameter of Q.2 inch and an exit diameter of 0.38 
inch. The locations of three points A, E, and X for comparisons are indicated here. 
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2.2 NOZZLE EXIT CONDITIONS 


Representative data on combustion chamber conditions and gas compositions were also obtained 
from Field (2001). The combustion gas is approximated as carbon dioxide (CO 2 ). Initial isen- 
tropic calculations, ignoring boundary layer effects, showed that a normal shock would be 
formed downstream of the nozzle exit. The nozzle exit conditions as derived from isentropic ex- 
pansion in the nozzle are shown in table 1. 

3. ANALYSIS 

3.1 CFD/KIRCHHOFF ANALYSIS 

The CFD/Kirchhoff analysis is based on the application of a CFD code for identifying the noise 
sources in the source-field and Kirchhoff surface integral for the propagation of sound radiation 
to the near-field and the far-field. In particular, the OVERFLOW CFD code (Buning et al„ 

2001) of NASA Langley Research Center (LaRC) is used for computing the instantaneous flow- 
field. This code solves the three-dimensional compressible turbulent flow Navier-Stokes equa- 
tions in generalized coordinates. A one-equation turbulence model of Spalart-AIlmaras (1992) is 
considered. The OVERFLOW code, in its present form, is designed primarily for the prediction 
of steady or unsteady turbulent flowfield and is used widely to model aerodynamic flows. It does 
not however provide the necessary time-dependent boundary conditions for handling the reflec- 
tion-free acoustic propagation. Therefore appropriate modifications to the code were made to 
provide a set of time-dependent reflection-free boundary conditions (includes periodic inflow, 
outflow, and radiation). 

For the acoustic radiation, the Kirchhoff code YORICK developed by Pilon and Lyrintzis (1998) 
of Purdue University was considered. The Kirchhoff surface, enclosing the nonlinear source re- 
gion, is chosen in a region where the linear wave equation is valid. The Kirchhoff method does 
not suffer dissipation and dispersion errors that could be encountered when the near-field and far- 
field sounds are directly calculated with the CFD code. By restricting the use of CFD methods to 
the near-field region for source identification, computational requirements are greatly reduced. 

3.1.1 AXISYMMETRIC GRID SYSTEM. For the CFD computations, an axisymmetric grid of 
200x100 (200 grid points in the axial direction and 100 grid points in the radial direction) is con- 
sidered. Figure 2 shows a schematic of the computational grid. A grid length of 60 jet radii and 
a grid radius of 10 radii are considered. The grid is clustered in the axial direction near the noz- 
zle exit for handling shocks and other high gradients in the flow and also near the outer radius of 
the nozzle for resolving shear layers. 

3.1.2 CFD SOLUTION. Appropriate time-dependent boundary conditions were applied to en- 
sure reflection-free boundaries. An outflow boundary condition of Thompson type (Thompson 
1990) is applied at the outflow boundary, which maintains the mean static pressure at the ambi- 
ent value. An acoustic radiation boundary condition of Tam and Webb (1993) is applied at the 
lateral boundary. At the inlet, the flow variables are specified based on the nozzle exit condi- 
tions. The CFD solution converged after about 15,000 time-step iterations (the code is run in a 
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time-accurate manner) before a periodic state is established. The solution is obtained on an IRIX 
workstation (SGI Indigo machine). About 30 hours of computing time are required for the solu- 
tion to achieve a periodic state. 

3.1.3 ACOUSTIC SOLUTION. After a periodic state is established, the appropriate data from 
CFD is communicated to the Kirchhoff code. The Kirchhoff surface is a cylindrical surface co- 
inciding with a grid line. The radius of the Kirchhoff surface is taken as about 6 radii from the 
jet axis. Only the lateral surface of the cylinder is taken into account ignoring cylinder ends due 
to the effects of nonlinearities. Specifically, the data to be specified on the Kirchhoff surface in- 
clude the instantaneous pressure, the pressure-time derivative, and the pressure gradient normal 
to the surface. The Kirchhoff code then computes the pressure-time signals and the OASPL in 
the near- and the far-field. From the pressure-time signals, it is possible to compute the SPL 
spectrum at any location with the aid of Fast Fourier Transform (FFT). Inside the Kirchhoff sur- 
face, the OASPL and the SPL spectrums are computed directly from the CFD solution for the 
instantaneous pressure. A more detailed account of the CFD and the Kirchhoff methods is pro- 
vided in a report by Kandula and Caimi (2001). 

3.2 CORRELATION/ANALYTICAL METHOD 

Leneman (1973) proposed an empirical correlation for the OASPL as 

OASPL (dB) = 115 + 10 log(FI) - 20 log(x) (la) 

where F denotes the total thrust (lbf), I the specific impulse (sec), and x the axial distance (ft) 
in the far-field from the jet exit plane. Margasahayam and Caimi (1999 and 2000) developed an 
analytical model to predict the far-field noise generated by rockets in the holddown state, repre- 
sentative of static tests. The OASPL is expressed by an extension of Leneman equation as 

OASPL (dB) = 115 + 10 log(FI) - 201og(x) - 20xl0 (x/50) (lb) 

The fourth term on the right-hand side is the near-field correction term due to Margasahayam and 
Caimi (1999). This equation represents a first-order estimate only and does not account for the 
effects of spectral composition, directivity, etc. 

The NASA SP-8072 method (1971) is considered in the rocket noise program (RNP) of Mar- 
gasahayam and Caimi (1999), which represents a more accurate correlation compared to the Le- 
neman correlation. It computes OASPL and SPL spectrums by an empirical correlation based on 
the normalized relative power spectrum as a function of Strouhal number derived from a number 
of solid- and liquid- fueled chemical rockets. 

3.3 PRESSURE FLUCTUATION SCALING METHOD 

The method of pressure fluctuation scaling due to Norton et al. (1999), which has been success- 
fully used for industrial gas piping networks, industrial gas turbine exhausts, blowdown systems 
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in liquefied natural gas (LNG) plants, and offshore platforms, etc., is discussed in the following 
paragraphs. 

The fluctuating pressure spectrum at the nozzle exit (after the shock wave) is estimated via a 
pressure fluctuation scaling technique based on (1) laboratory data for a turbulent wall pressure 
spectrum in a cylindrical shell, (2) a turbulent wall pressure spectrum on a flat plate, and (3) a 
fluctuating wall pressure spectrum downstream of a severe internal flow disturbance such as a 
depressurization valve or a 90-degree mitre bend and using the post-shock centerline input data 
(static pressure, static temperature, density, flow velocity, and Mach number) from the CFD re- 
sults. The nondimensional fluctuating pressure spectral density <J> P is given by 


®,(Q) 


4 G p m 

p 2 U 3 a 


where the nondimensional frequency £2 (Strouhal number) 


Q = (DalU 


(2a) 


(2b) 


where G p is the fluctuating pressure spectral density, p the gas density, U the flow velocity, a the 
nozzle exit radius, and co = 2 nf . Controlled laboratory measurements of G p , p, and U for a 
given turbulent boundary layer flow in a shell of internal radius a are used to calculate the non- 
dimensional fluctuating pressure spectrum <J> P according to equation (2a). The frequency axis is 
simultaneously scaled according to equation (2b). Hence, a spectrum of nondimensional fluctu- 
ating turbulent pressure spectral density d> p against nondimensional frequency Q is obtained. 

The turbulent wall pressure spectrum (shell or flat plate) provides a lower level, and the severe 
internal flow disturbance provides the upper level for scaling purposes. 


The sound spectrum is then calculated at the exit downstream of the shock. Power spectral esti- 
mates of fluctuating pressure against radian frequency for the tabletop rocket are obtained from 


S,W 


® p (Q)p 2 U 3 aAco i 

4 


(3a) 


where co i =UQ.,la (3b) 

and A(Oj is the nominal bandwidth of band i with center frequency (Oj. 

Correction factors are then applied for the open jet exhausting out of the nozzle. At this stage, a 
1/r decay with distance from the nozzle is assumed for circumferential radiation away from the 
exit plane (e.g., position A, figure 1). The 1/r decay is also considered for post shock centerline 
decay (e.g., position X) because the mean flow will enhance/convect the monopole- and dipole- 
type radiation from the duct exit. These assumptions need to be refined but serve as an order of 
magnitude prediction at the present time. 
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The results were computed in linear octave frequency bands and in spectral values [decibel (dB) 
linear/hertz (Hz)]. Both the lower- and the upper-limit scaling spectra were utilized to estimate 
the pressure fluctuations downstream of the shock wave at all three locations downstream of the 
nozzle exit. The pressure fluctuation scaling results presented here thus provide both lower- and 
upper-level estimates at locations A, E, and X for completeness. 

4. RESULTS AND COMPARISON 

4.1 FLOWFIELD RESULTS 

The axial velocity contours and the temperature contours from CFD are displayed in figures 3 
and 4 respectively. An overview of the jet flowfield structure, including the shear layer mixing, 
is evident from these plots. A normal shock is formed near the exit plane, thus resulting in a sub- 
sonic flow downstream of the exit. The existence of a normal shock outside the nozzle is to be 
expected based on the isentropic expansion/normal-shock calculations for the nozzle with the 
given area ratio (Mach number) and chamber to ambient pressure ratio (Shapiro 1953). The 
temperature contours suggest the formation of vortical structures in the shear layer. 

Figures 5 through 8 represent the centerline variation of important flow parameters (static pres- 
sure, static temperature, flow velocity, and Mach number). Figure 5 displays typical static pres- 
sure excursions in the core region and the pressure approaches the ambient value far downstream 
of the nozzle exit plane. It is inferred from figure 7 for the centerline axial velocity distribution 
that the core length of the jet extends to about 25 jet radii from the exit plane. As evident from 
figure 8, the average Mach number in the core region downstream of the shock is about 0.6, 
characterizing the jet as subsonic. 

4.2 ACOUSTICAL RESULTS 

Figure 9 presents the OASPL contours in the near- and the far-field {xlr j from -60 to 60, and 

z/ r } from 10 to 70) as computed by the CFD/Kirchhoff formulation. The OASPL ranges from 

80 to 130 dB in this region. Two primary lobes are recognized and located downstream of the 
exit plane. The absence of a dominant Mach wave radiation is due to the subsonic nature of the 
jet. The OASPL contours within the jet, characterizing primarily the nonlinear source-field, dis- 
played in figure 10, are evaluated directly from the CFD solution. Inside the jet shear layer, the 
sound pressure levels are high (as much as 180 dB), with the highest OASPL occurring about 7 
jet radii downstream of the exit. The relatively high-sound levels at this position are attributable 
to the high-turbulence intensities in the shear layer. Figure 11 displays the OASPL distribution 
along the line BC, 24.52 radii from the jet axis (figure 1), as computed from CFD. The OASPL 
is seen to be highest near the exit plane and decreases away from the exit. 

In the following plots, comparisons are made for the OASPL and spectral distribution of SPL at 
three representative positions A, E, and X (see figure 1). Point A (jt/r y =0, and z/r =24.5) is 

located in the exit plane, 24.5 jet radii away from the jet axis; point E {x!r j =0.63, and z!r j =0) 
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is situated on the jet axis, 0.63 jet radii downstream of the exit plane; and point X (jc//v=23.5, 

and z / Tj =1) is stationed 23.5 jet radii downstream of the exit plane, and one jet radius away 

from the axis. Table 2 summarizes the comparison of OASPL at stations A, E, and X as esti- 
mated by the different methods. 

The spectral composition of 1/3 octave SPL at point A in the exit plane is presented in figure 12. 
A similarity is noted in the SPL spectrum from the CFD and the upper-level PFS, although the 
SPL from the upper-level PFS exceeds the CFD result. However, both the CFD and the upper- 
level PFS method suggest a peak frequency of about 20 kilohertz (kHz). This peak frequency 
corresponds to a Strouhal number St=0.20. The lower-level PFS method produces SPL that is 
close to that of the CFD but produces a peak frequency of about 100 kHz, which appears some- 
what overestimated. The deviation in the SPL between the upper-level and the lower-level PFS 
is attributable to the fact that the lower-level result is based on turbulent pressure spectrum for 
shell/flat plate, whereas the upper-level result is based on the pressure spectrum scaling corre- 
sponding to severe internal flow disturbance, as indicated before. High value for the peak fre- 
quency level (of the order of kHz) for the SPL distribution is a characteristic of the small-scale 
nozzle considered here, as to be expected from the similarity scaling provided by the Strouhal 
number. 

A comparison of the 1/3 octave SPL content at location E on the axis downstream of the exit 
plane is portrayed in figure 13. In general the SPL values at E are considerably higher than those 
at position A, due to the near-field location of E and the effect of directivity (E is located in the 
jet and is also closer to the exit plane, whereas A is located in the exit plane and oriented at 90 
degrees to the jet axis). A peak frequency of about 35 kHz (corresponding to St=0.35) is indi- 
cated by CFD, the NASA-SP method, and the PFS upper-level method. The CFD results for the 
SPL are higher than those from the NASA-SP method. In the case of the PFS method, the upper- 
level values for the SPL are closer to the CFD, whereas the lower-level results match well with 
the NASA-SP method. 

Figure 14 compares the spectral distribution of 1/3 octave SPL at location X, which is down- 
stream of location E but one jet radius away from the axis. In general the sound pressure levels 
at X are smaller than those at E. A reduced peak frequency of about 25 kHz (corresponding to 
St=0.25), relative to location E, is indicated by CFD, the NASA-SP method, and the PFS upper- 
level method. The SPL variation between the various methods at position X is somewhat smaller 
as compared to that at positions A and E. 

A comparison of the 1/3 octave SPL at locations A, E, and X as predicted by CFD is shown in 
figure 15. Downstream of the exit plane, the peak frequency level is increased in the near-field, 
as the distance from the exit plane is decreased. While there appears to be a similarity in the 
shape of the SPL spectrum between E and X, there is a departure in similarity at position A for 
frequencies in excess of about 40 kHz, where a rapid decline in SPL is observed. This lack of 
similarity in the SPL spectrum at position A is perhaps attributable to the existence of a normal 
shock near the nozzle exit or to the far-field directivity effects. 
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5. CONCLUSIONS 


The CFD/Kirchhoff method provided flowfield structure and distributions of the OASPL and 
SPL spectrum for the tabletop rocket. As a result of the overexpansion in the nozzle, a normal 
shock forms in the immediate downstream of the nozzle exit, providing a subsonic jet. Vortical 
structures in the jet shear layer are observed. High values of OASPL (140 to 180 dB) are noted 
in the shear layer region comprising about 15 jet diameters downstream of the exit plane. Lower 
levels of OASPL (80 to 130 dB) are computed in the near- and the far- field region. Because of 
the very small size of the nozzle exit diameter, the jet is distinguished by high peak frequency 
levels (20 to 35 kHz) for the SPL spectrum. 

A comparison of the results from CFD, the NASA-SP method, and the PFS method at three spa- 
tial locations (A, E, and X) showed that a similarity in shape is generally noted for the SPL com- 
position. There is some disagreement between the different methods for the peak frequencies 
and the SPL values. In the jet region (positions E and X), the CFD predictions for the SPL are 
found to be higher than those from the NASA-SP method. The upper-level results in this region 
agree better with the CFD results, whereas the lower-level results match better with the NASA- 
SP values. In the far-field region (location A), the lower-level PFS compares better with the 
CFD results, while the upper-level PFS values considerably exceed the CFD predictions. A more 
realistic assessment of these methods can only be carried out with the aid of an experimental test 
program currently under development. 
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Table 1 . Summary of Nozzle Parameters 


Parameter 

Value 

Stagnation pressure, psia 

54.7 

Stagnation temperature, °R 

1960 

Nozzle mass flow rate, lbm/s 


Nozzle throat diameter, inch 

■ 

Nozzle exit diameter, inch 

0.38 

Exit pressure, psia 

2.38 

Exit temperature, °R 

987 

Exit velocity, ft/s 


Acoustic velocity at exit, ft/s 


Nozzle exit Mach number 

2.65 

Nozzle exit density, lbm/ft 3 

0.0099 

Molecular weight for CO 2 

44 

Specific heat ratio for CO 2 

1.28 

Jet exit Reynolds number 

2.9xl0 4 

Ambient pressure, psia 

14.7 

Ambient temperature, °R 

519 


Table 2. Comparison of OASPL at Locations A, E, and X 


Location 

CFD 

Leneman 

NASA- 

SP-8072 

PFS 

Original 

Modified 

Upper 

Level 

Lower 

Level 

A 

118 

_* ** 

- 

- 

153 

119 

E 

151 

175 

155 

125 

176 

130 

X 

140 

144 

124 

117” 

135 

122 


* Not applicable 

** Value corresponding to z = 0 
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Figure 2. Axisymmetric Grid (200x 100) for CFD 
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Figure 9. Contours of OASPL Predicted From CFD/Kirchhoff Method 
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Figure 12. Comparison of 1/3 Octave Spectra at Location A. 



Figure 13. Comparison of 1/3 Octave SPL Spectra at Location E 
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Figure 14. Comparison of 1/3 Octave SPL Spectra at Location X 



Figure 15. Comparison of 1/3 Octave SPL Spectra From CFD at A, E, and X 
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